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Abstract. The phases and excitation spectrum of an easy-axis ferromagnetic 
chain of S = 1/2 magnetic impurities built on the top of a clean metallic surface 
are studied. As a function of the (Kondo) coupling to the metallic surface and 
at low temperatures, the spin chain exhibits a quantum phase transition from 
an Ising ferromagnetic phase with long-range order to a paramagnetic phase 
where quantum fluctuations destroy the magnetic order. In the paramagnetic 
phase, the system consists of a chain of Kondo-singlets where the impurities 
are completely screened by the metallic host. In the ferromagnetic phase, the 
excitations above the Ising gap are damped magnons, with a flnite lifetime arising 
due to the coupling to the substrate. We discuss the experimental consequences 
of our results to spin-polarized electron energy loss spectroscopy (SPEELS), and 
we finally analyze possible extensions to spin chains with S > 1/2. 
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1. Introduction 

Recent progress in scanning tunneling microscopy (STM) and other surface- 
manipulation techniques at the atomic scale has made it possible to build low- 
dimensional magnetic structures with unprecedented control [TJ Atomic-scale 
magnetic structures possess an enormous potential for applications in spintronics as 
well as in classical and quantum computation due to their ability to store and process 
(quantum) information in the atomic spins. However, in order to build reliable atomic- 
size memories or qubits out of such systems, it is crucial to understand how the 
coupling to the environment affects their properties. For instance, it is known that 
quantum nanomagnets coupled to electronic degrees of freedom are prone to the effects 
of decoherence induced by Ohmic quantum dissipation |3) |4J . 

In addition to their interest for information technologies, low-dimensional 
magnetic systems coupled to dissipative environments are also very important from 
the point of view of fundamental physics. This is because they allow to study the 
crossover from quantum to classical behavior as the coupling to the environment 
is increased ^ |S] . Current theoretical models that describe local quantum degrees 
of freedom coupled to sources of dissipation predict rich and interesting physical 
phenomena at low temperatures, including quantum phase transitions and quantum 
criticality |BHO]j as well as broken continuous symmetries and long-range order (LRO) 
in low-dimensional systems |10HT5] . The interplay between quantum fluctuations, 
which are induced by quantum confinement and reduced dimensionality, and quantum 
dissipation, caused by the interaction with the environment, is at the core of this 
exotic behavior. From this perspective, magnetic nanostructures built on the top of 
clean metals offer a unique platform to test our understanding on dissipative quantum 
systems. For instance, the simplest zero-dimensional (OD) case of a single magnetic 
impurity coupled to a metallic environment [such like a Co or Mn adatom sitting 
on a Au(lll) surface] is a practical realization of the celebrated Kondo effect, one 
of the most paradigmatic phenomena in the physics of strongly correlated electrons 
[16] . This phenomenon is the spin-compensation by the conduction electrons of the 
magnetic moment of an impurity embedded in a metallic host at temperatures 
T ^ Tk, where Tk is the so-called Kondo temperature. Eventually, at T = 0, the 
ground state of the system is a many-body singlet (termed a 'Kondo singlet'), formed 
by a linear superposition of and the spin of a cloud of conduction electrons that 
is anti-ferromagnetically correlated with the impurity local moment. The original 
magnetic impurity is thus completely screened by the metallic environment. In 
the context of magnetic adatoms deposited on clean metallic surfaces, this is a well- 
established phenomenon, which reveals itself as a narrow Fano resonance at the Fermi 
energy in the STS spectra \^\^^ . 

Unfortunately, the extension of the single Kondo impurity to the case of many 
interacting magnetic impurities embedded in a metallic host represents a formidable 
theoretical challenge. Due to the many-body nature of the problem, exact methods 
are only limited to simple cases (i.e., consisting of few magnetic impurities), and at 
the cost expensive numerical calculations [20'. On the other hand, dynamical mean- 
field theory (DMFT) methods [21, provide reliable results in the case of bulk 3D 
compounds, but are much less successful when applied to low-dimensional systems. 

One particular class of such low-dimensional systems, namely, one-dimensional 
(ID) chains of magnetic atoms or organic molecules with active magnetic atoms, are 
currently under intensive experimental investigation |22H26j . Thus, by depositing 
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~ 0.1 monolayer (ML) of Co on Pt(997), which is a vicinal surface containing 
steps, chains of magnetic atoms have been recently created [22]. The magnetic Co 
atoms decorate the steps and form ferromagnetic chains, which, when probed by an 
external magnetic field, typically exhibit superparamagnetism. However, in the case 
of Co chains at temperatures T < 8 K, the magneto-crystalline anisotropy results 
in relaxation times which are longer than the experimental typical timescale, and 
the response of the system becomes effectively ferromagnetic |22| . Another recently 
reported system, consisting of chains of Co atoms on Cu(775), has been shown to 
undergo a Peierls distortion leading to dimerization |57]. This dimerized phase has 
been explained by the strong correlations existing in the partially filled d shells of 
the Co atoms, which are in the high spin configuration. Maintaining the high spin 
in a geometry where the overlap between the d orbitals is small and correlations are 
strong, favors a ferromagnetic ground state, which in turn makes the charge-density 
wave instability possible. 

In the above-mentioned systems, the coupling to the metallic substrate was 
assumed to be weak, and therefore was neglected in their analysis. However, as 
we show below, the metallic environment can have dramatic effects for the phase 
diagram at low temperatures, and it cannot be neglected when studying the magnetic 
excitations of the chain. With these systems in mind, in this contribution we focus on 
the analysis of an easy-axis ferromagnetic (FM) ID spin chain deposited on the top of 
a metallic surface. Interestingly, as a function of the coupling with the substrate, we 
predict a quantum phase transition from an Ising ferromagnetic phase with long range 
order to a paramagnetic phase where quantum fluctuations proliferate and destroy the 
ordered phase. Although modifying the coupling to the substrate in the experimental 
systems created so far may be hard, in the future other systems can be created 
where the coupling to the substrate could be tuned, uncovering even more interesting 
phenomena. 

The article is organized as follows: In the next Section, we introduce a simple 
model to describe an easy-axis S" = 1/2 ferromagnetic chain on a metallic substrate. 
This model is studied in the framework of mean-field theory in Section [3] and the 
quantum phase diagram is thus obtained. The phase diagram is displayed in Fig. [4] 
and it is one of the most important results of this work. In Section |4] we consider the 
excitation spectra of the ferromagnetic phase. We show that a weak coupling with 
the metallic substrate can have important consequences for the magnetic excitations 
of the ferromagnetic chain. This properties of the excitations, and in particular, the 
substrate-induced damping, can be accessible by experimental probes such as spin- 
polarized electron energy loss spectroscopy (SPEELS). In Fig. [5] we plot the spin- 
response of the ID spin chain due the substrate-modified magnon excitation spectrum, 
which is also an important result of this work. Finally, in Section |6] we offer our 
conclusions and outlook for future research directions. In the [Appendix A[ we give 
more details about the technical aspects of our calculations. 

2. Model 

Let us consider a chain of magnetic atoms deposited on top of metallic surface (see 
Fig. fl]). We model the magnetic interactions between the atoms are using the spin-^ 
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FM spin chain 




Figure 1. Graphical representation of a FM spin chain deposited on the top of 
a metallic surface. Kondo-screening of the individual spins at low temperatures 
T Tfc tends to destroy the FM long-range order and induce a quantum phase 
transition to a paramagnetic phase (cf. Fig. Hjl. 



XXZ Hamiltonian: 

N 

-ffo = - ^ [Jh {S^Sf_^i + S'f_^iS'f ) + JHS^Si_f_i] . (1) 

i=l 

In the above expression, we have assumed anisotropic exchange interactions, which 
originate from the reduced symmetry of the surface crystal environment, as well as 
from spin-orbit interactions of the electrons in the metallic substrate [2]. Assuming 
that the impurity spin is = ^ implies that the single-ion anisotropy DJ^ii'^i)^ 
reduces to a constant term. Furthermore, for S' = exchange terms like 
(Si • Si+i)^ can be expressed in terms of (Si • Si+i). However, for impurities 
with S > such terms must be taken into account, but we shall not pursue their 
analysis here. 

The nature of the ground state of Eq. (jlj depends crucially on the ratio A ~ 
Jfj I Jj^ , where we assume jjj^Q and Jfj > 0. For A > 1 (A < -1), Eq. ^ describes 
an easy-axis (i.e. Ising-like) ferromagnet (anti-ferromagnet) . In this system, the spin 
waves (magnons) and their bound states are separated from the ground state by a 
gap p8l I29|. On the other hand, for |A| < 1, the spin chain Hamiltonian ([ij exhibits 
a XY phase with a gapless spin- wave spectrum [30, 31 . However, as anticipated in the 
introduction, in what follows we will focus only on the FM Ising limit A > 1 of Q, and 
we refer the reader to Ref. |15j for a complementary study in the case of easy-plane 
anisotropy, where the spin chain exhibits an XY phase. Nevertheless, we note that 
the results reported in Section |3] are also applicable to the antiferromagnetic (AFM) 
regime corresponding to A < —1 (see Section [6]). Finally, for the sake of simplicity, 
we have neglected longer-range magnetic exchange interactions in Eq. ([ij, may be 
present due to the Rudermann-Kittel-Kasuya-Yosida (RKKY) mechanism mediated 
by the metallic substrate [35;. We note that these interactions can be easily included, 
as long as they do not modify the nature of the Ising FM ground state. 

The coupling between the spin chain and the metallic substrate is described by 
an anisotropic Kondo-exchange Hamiltonian |16] : 

N 

Hk =Y.Jk [s'ts'' (RO + sfsy (R.)] + rkSls^ (Ri) , (2) 
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Figure 2. Local bath approximation. For impurities embedded in a 2D or 
3D metallic substrate, the interference of Kondo screening clouds is negligible 
when the impurities are separated by a few Fermi wave numbers kp^ |36H38| , and 
therefore can be considered as independently screened (ISj. 



where s° (Ri) (a — x,y,z) is the spm-density of conduction electrons at the position 
of the z-th spin impurity in the surface Rj — i.ao (with oq the lattice parameter of 
the chain). Note that this model bears some similarity with the Kondo-lattice model, 
which is used to describe the properties of heavy- fermion systems |T6l [33| [34] . In 
the case of impurities with S" > | , the Kondo exchange involves more channels (with 
channel-dependent couplings) and results in a dynamical screening of the impurity 
spins occurring in several stages, involving different Kondo temperatures |35 . We 
shall not pursue such an analysis here. 

Despite its simplicity, the model introduced above still represents a formidable 
theoretical challenge, and in order to make analytical progress we shall introduce 
some simplifying assumptions. We approximate the metallic substrate by a collection 
of one-dimensional metals (cf. Fig. [2]) described by the Hamiltonian 

N 

where c^^ (fc) creates an electron with linear momentum k and spin projection a. The 
index i indicates that the electron reservoir described by the operators cj^.^ (fc) , Cai (k) 
is only coupled to the impurity located at Ri and not to others with j ^ i- For this 
reason, we dubbed this approximation the "local bath approximation" (LBA) |15| . 
Within the LBA, the spin-density operator reads: 

s'^(R,) = 5:cL(fc)f^) c,,(fc), (4) 

a,/3 \ ^ / 

where cr" are the Pauli matrices. 

The dimensionality of the substrate is crucial for the justification of the LBA. 
When the magnetic impurities embedded in a high-dimensional (2D or 3D) metallic 
substrate are separated by a distance of the order of a few k^^ , with kp the Fermi 
wavevector of conduction electrons, the interference of Kondo screening clouds that 
dynamically quench the magnetic moments, becomes negligible |36fl38| . Consequently, 
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at distances oq ^ kp^ the impurities can be considered as independently screened 
|15| . This is an important difference with respect to the case of a strictly ID metallic 
substrate, as in the case of the ID Kondo-lattice model [32] ■ In that case, the single 
Kondo-impurity limit is reached only at distances ao ^ S,k vf/Tk- Since the 
Kondo temperature Tk is an exponentially small scale, in the purely ID geometry the 
single-impurity regime is only reached at extremely dilute impurity-spin concentrations 
[36H38]. 

In actual experiments, although the magnetic nanostructures are built on the top 
of a metallic surface, kp is determined by the 3D bulk electron density due to a non- 
vanishing overlap between the bulk and surface conduction states [15]. For instance, 
the Fermi wave number for Au bulk free conduction electrons is 7r/fc^""^ — 2.6A, while 
for Au(lll) surface states is tt/A:™''^ = 16. 6A [40]. This fact dramatically increases the 
range of applicability of the LBA. The physical picture provided by the LBA is also 
supported by the behavior of the STS Fano line shapes in experiments on magnetic Co 
atoms deposited on Cu(lOO) and separated by distances uq > 8 A, which are identical 
to the single-impurity STS line shapes |41 j . 

In addition, although in the actual experimental realizations the limit oq ^ kp^ 
might not be strictly realized, the model resulting from applying the LBA to Eqs. ( T|2 1 



can be considered as the simplest "toy model" that captures the competition between 
the Kondo physics, which is a local quantum critical phenomenon, and the magnetic 
interactions along the chain [TS]. Concerning the latter, it is worth noting that the 
LBA cannot reproduce the non-local RKKY exchange coupling mediated by the metal 
and arising at order O {Jx:)- As mentioned before, this is not a crucial drawback, 
since one can always redefine the couplings J^tJh ^'i- ^ske the RKKY 
contribution explicitly into account. More importantly, what eventually justifies the 
separate treatment of RKKY and Kondo interactions is that, although they both 
originate in the same Hamiltonian the RKKY interaction results from electronic 
states deep inside the Fermi sea and is a static coupling, while the Kondo effect is 
purely dynamical effect originated predominantly in the Fermi surface |34l 112] . 

3. Mean-field analysis of the T — phase diagram 

In order to gain insight into the low-temperature phase diagram of the system, we 
now focus on the limit A > 1 of Eq. ([T]), describing an isolated FM Ising spin 
chain. In that case, the ground state of the isolated chain is separated from the 
single-magnon excitation spectrum by a gap Eg = Jfj (l — A~^) and by a gap 
Eg = Jfj (^1 — A~^) > Eg from the two-magnon bound states [^ This means 
that the long-range FM order in the isolated chain is stable at low temperatures and, 
to a first approximation, it is safe to take = in Eq. ([I]). 

We now introduce the coupling to the substrate Eq. ([2]), which has important 
consequences for the low-temperature behavior of the the spin chain. For instance, the 
nature of the ferro magnetic to pa ramagnetic phase transition is profoundly modified. 
As we show in the Appendix A in the case — the model for the Ising chain 
coupled to the substrate can be mapped onto the ID dissipative quantum Ising model, 
which is in the universality class of the classical 3D (instead of ID) Ising model. In 
the following, we exploit this fact to introduce a mean-field (MF) approach, which 
therefore provides a good description of the phases of the spin chain in the (easy-axis) 
Ising limit. 

We introduce the MF approximation to the Ising term in Eq. ([T|) by making the 
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replacement SfSf^^ ^ Sf (5f+i) + (S'f ) S'f+i - (S'f ) (5f+i), which leads to the MF 
Hamiltonian of the spin chain: 



Hr-'^-T.K.s:, (5) 

where 

/ieff - 2J^ (Sf ) (6) 

is the local Weiss field. The MF Hamiltonian of the system is therefore written as 

HMF^Y.nr + ^. (7) 

nf"^ = [S'ts- (R.) + Sfsy (R,)] + JliSts'- (R,) 

— hcffS- + T-Lpi, (8) 

with Hpi = "^ka^k c]ji (fc) Ccti (fc) describing the i-th conduction electron bath. In 
the following we drop the lattice sub-index i in Hf^^ and Hpi in order to lighten the 
notation. 

We note that the MF approximation together with the LBA leads to a set 
of independent Kondo-impurity problems in an effective magnetic field hcff, which 
is a quantity that must be obtained self-consistently. The MF thermodynamical 
properties of the spin-chain can therefore be obtained from the free-energy / (T, h^s) = 



-TlnTr 



for a single Kondo impurity in a magnetic field /icff- Note 
that the Weiss field hcs does not act upon the conduction electrons, it only acts on 
the impurity. The mean-field self-consistency condition is obtained from the relation 
between the impurity magnetization, 

JTlimp \T, /leff) = {S ) = — (9) 

and the definition of the Weiss field Eq. ([6|: 

minip (T,/loff) = (10) 

To make further progress we now need to know the function mimp (r, h^s) for a 
single Kondo impurity. Introducing the Abelian bosonization method |3T] to describe 
the metallic ID chains in Fig. [2j and the unitary Emery-Kivelson transformation 
143] [cf. Eq. ( A.3|] , the anisotropic Kondo model can be mapped onto the resonant- 



level model |44H46] . The resonant-level model is exactly solvable when the parameter 
is chosen at the so-called "Toulouse line" |37], defined by tan~^ {T^Pobo^K f"^) — 
I ^1— '^here po and bo are the density of states and lattice parameter of 



the bosonized ID chains, respectively (cf. Appendix A I. In combination with the 
renormalization group (RG) approach [48 , which describes the fiow of parameters 
{ (^) 7 J^K (^) } of anisotropic Kondo model as a function of the conduction- 
band cutoff parametrized by A {£) = AqC"^, the mapping to the resonant-level model 
allows to describe the physical properties of the strong coupling limit of the Kondo 
model. One starts with the "bare" couplings { J^, J^} = {J^ (^ = 0) j Jr = 0)}, 
and renormalizes up to the point where (i* ) — , and at that point we use the 
exact solution of the resonant-level model jST]. However, we stress that mimp (T, /icff) 
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Figure 3. Graphical solution of the mean-field equation, Eg. |10| for T = 0. For 
Jfj larger than the critical value there are two kinds of mean-field solutions: those 
for which hefl < ^^nd those for which h^ff > Tj^. In the first case, the coupling 
between the impurities does not suffice to destroy the Kondo singlet. However, in 
the latter case, the coupling is strong enough to break the Kondo singlets. 



can be obtained in a more general case (and with more sophistication) using any 
"impurity-solver" method, such as the numerical renormalization group (NRG) |16] . 
The mapping to the resonant-level model yields the result |44fl46] 



QP {T, /icff) = -Im 

TT 



K 



ihcs 



(11) 



2ttT 2nT, 

where ip (z) is the digamma function ^49) and the broadening of the resonant level 
r = pobo I 4^ 1 can be interpreted as Tk ^ F. 

The phase boundary between the FM phase where (S^) cx h^e ^ 0, and the 
paramagnetic phase with (S*^) = is obtained from (10 1 in the limit /left (cf. 
Fig |3]). For \T iTK^h^slTK] ^ 1, the impurity exhibits a characteristic (Pauli) 
paramagnetic response mimp (T, /leg) = Ximp (?") ^leffj where 



Ximp (r) 



1 ,y/l Tk 



Ximp (0) 



1 - 



T 
Tk 



O 



TV 
Tk 



(12) 
(13) 



where Ximp (0) = I/t^Tk is the magnetic susceptibility of the impurity at T = 0. On 
the other hand, for /icff S> min {Tr-, T}, TOi„ip (T, hcfi) approaches the saturation limit 
of 1/2 as 



rriimp (T, /icff) — — arctan 

TT 



27rT + Tk 



(14) 



The above result can be also obtained using the scaling equations, and it is modified 
by anisotropy at the Toulouse point. Thus, a non- vanishing solution for /loff is found 
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Figure 4. Mean-field phase diagram of the Ising chain coupled to the metallic 
substrate. At low enough temperatures, and in the regime Jf^ <g; Ti^ quantum 
fluctuations induced by the (Kondo) coupling to the substrate completely destroy 
the long-range magnetic order. 



from Eq. ( 10 1 in the limit hcff 
1 



2J 



< Ximp (0) 



H 



1 



11 ( 2L 



(15) 



From the above considerations, the phase boundary is therefore given by the critical 
line 



Tk 



V3 



2J 



(16) 



H 



and from here we can obtain the MF phase diagram of the spin chain (See Fig. |4|. 

It is interesting to note that, in the FM phase, as Jfj/TK increases from the 
phase boundary, the value of the Weiss field /leff crosses over from a regime /lofi <C 
where the Kondo singlet is still robust, to the regime of large hes ^ Tk where the 
Kondo singlet is destroyed (cf. Fig. |3|. Therefore, the above mean-field theory yields 
a regime where the Kondo effect and a FM LRO coexist. In the following Section, we 
shall return to this point again. 



It is also worth mentioning that the result in Eq. (16) clearly contrasts with the 
case of the isolated Ising chain, for which the analytical form of Tc Eq. (16 1, is very 



different from the expression for Tc encountered in the study of the classical Ising 
model, where Tc oc Jfj [50 . Physically, in the present case, the spin-flip processes 
induced by the Kondo coupling to the metal Eq. ([2| can be considered as quantum 
fluctuations that tend to destroy the LRO in the Ising chain as T — ^ for Jfj/Tx ^1-6 
(See Fig. |4]). 

In the following Section we characterize the excitations of the ordered FM phase. 
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4. Magnetic excitations in the FM phase 



We next focus on the FM phase, in the regime J^/Tk <^ 1 (right-bottom region in 
Fig. |4| , where the spectrum of the system is gapped and the effective Weiss field hes 
is large enough so as to inhibit the Kondo screening of the impurities. We shall focus 
in particular on single magnon excitations that can be created by probes (such like 
spin-polarized electron energy loss spectroscopy, SPEELS, see Section |5| that couple 
weakly to the ferromagnetic chains. Such probes are assumed to produce single spin- 
flips, and thus the calculation of the system response can be carried out within linear 
response theory. This approximation allows us to neglect the effect of magnon bound 
states in the spectrum of the FM XXZ chain. Thus, the low-lying excitations involving 
a single spin-flip can be described using the Holstein-Primakoff (HP) representation 
of the spin operators [51] : 



S+ ::^b]\ 2S-b]bi = V2S 



2S 



b]hb. 



'2S 



b\ + 










6. + ( 







Sf = b]b, - s 



(17) 

(18) 
(19) 



where bi,bl are bosonic annihilation and destruction operators, obeying 
which physically represent a single spin-flip. For the XXZ chain introduced in Sec. [2] 
we should take S — ^in the above expressions. However, below we shall work assuming 
that S is large, which is known to provide an accurate description of the lowest energy 
excitations [51 . Assuming S large means that our results also apply to spin chains 
with S* > ^, provided the Hamiltonian is of the XXZ form and we assume that 
the Kondo-exchange coupling with the substrate describes the strongest channel of 
substrate electrons that couple to each impurity along the chain. Thus, working to 
leading order in S (i.e. in the large-S* approximation), we find: 

N 

J2[-Jhs'-Jks{blh+i + bl, 



HP 



{blb,+4+,b,, 

which can be easily diagonalized in a running-wave basis 



(20) 



(21) 



and leads to 



= E,b% - NJ}jS' + O (1) , (22) 

where 



9 



Eg ^ 2SJfj - 2Sjjj cos (qao) , 

^ 2S (J^ - Jh) + 25 + O , (23) 

is the dispersion relation for magnons, and oq is the lattice parameter. In addition, 
in Eq. (22) we have neglected terms of O (1), which involve interactions between the 
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magnons. From Eq. (23 1, we see that for Jfj > > (i.e. A > 1) the excitation 



spectrum of the spin chain is gapped, and the spin waves exhibit a quadratic dispersion 
for momentum g — > 0. 

Let us next consider on the Kondo-exchange coupHng to the metal. Eq. We 
shall treat this coupling using the HP representation in the large-S* limit: 



- E [JK\/Sj^{bls- (R,) + s+ (R,) 6, 
+ (bib, - (R,) 



(24) 



As we will see, although this Hamiltonian no longer describe the Kondo effect (we recall 
that we are working in the limit of large h^s) , the coupling to the gapless degrees of 
freedom in the metal still can have important consequences to the spin chain. To 
make further progress, we now need to integrate out the electronic degrees of freedom 
in the conduction band. To this end, we shall rely upon the coherent-state functional 
integral representation of the partition function |52| : 



r ^ 



where [bi, 6^] is the effective Euclidean action of the spin chain 

sS,^[b,b] ^So"^[b,b]+s^,rsM, 



where the first term is the Euclidean action corresponding to ( 22 1 

[b. b] =J2 d^'^i (^) + ^^/fi] b, (r) , 
q Jo 

= ^ ^ [-iflUJ,i + Eg] bq (Wn) bq (w„) , 



(25) 

(26) 

(27) 
(28) 



and S^^^ [b, b] = — In / exp — L dr H^^ (r) \ is the contribution arising from 

\ L J / Fi 



Hamiltonian ( |24p , where we integrate over the fermions {cai,Cai). Note that this can 
be done exactly, as Hp is quadratic (we assume that Fermi liquid theory applies) and 
leads to 

c ^ [-hp 



i=l 



X [6,(r')MT')-^]}, (29) 
where we have introduced the spin correlation functions in the conduction band: 



[s+ (R„r)s- (R„0))f, 
[s' (R„t) (R„0))f. 



(30) 
(31) 
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Note that these correlation functions do not depend on because all the local ferniion 
baths are identical. Now. since these baths are assumed to be Fermi liquids, we have 
that at T = and for r ^ — (where ep is the Fermi energy of the metallic substrate) : 



A± 

r2 ■ 



Xfir)--'-^, XFW--^, (32) 

which follows from the existence of a ^ w (where uj is the excitation energy) spectrum 
of particle-hole excitations in the Fermi sea at low energies. Note that since Xf ^'^^ XF 
have units of {E x t)~^ (where E and t denote energy and time units, respectively), 
which stem from the prefactor of h in Eq. 30 and 31 then A* , A^^ have units t/E . In 

(33) 
(34) 



Fourier representation, we express (32 1 as: 
Xf i^n) ^ A=^ \uJn \ , 

where w„ = 27rn//3 are the Matsubara frequencies, and therefore Eq. (29) becomes 



N 

■EE 

w„5^0 i=l 



i^7i\bi {U!n)hi (w„) , 



(35) 



where we have dropped the last term in Eq. ( 29 1 because it amounts to a (retarded) 
magnon-magnon interaction. This is consistent with neglecting magnon interactions 



in (23 1. The above contribution is quadratic, and replacing it into (26) leads to the 



effective action 



qHP 



where we have defined the magnon propagator 

1 



9q i^n) 



(36) 



(37) 



with a cx A* (J-^ 



i^j a dimensionless parameter quantifying the dissipation caused by 
the conduction electrons on the spin chain. We next perform the analytic continuation 
to real frequencies, in order to obtain the retarded magnon propagator. Thus, 
replacing iw„ — > w + zO"*" and |aj„| — >■ —iuj, yields 

<H= .,. J , (38) 



huj — Eg + iauj ' 

Thus, the pole of the propagator becomes hwg = Eg {1 + ia/h)' 



E„ 



for 



small a, which implies that the coupling to the metal Eq. ( 24 1 induces a finite lifetime 
Ha~^ and damping on the magnon excitations. Physically, this happens due 



r 



to the presence of gapless particle-hole excitations in the excitation spectrum of the 
metallic substrate, and it can be interpreted as Landau damping [53-'--55] . Interestingly, 
this damping does not vanish in the limit g — > 0. At first sight, this may seem 
surprising, since in the g — > limit, magnon excitations correspond to a uniform 
magnetization along the z direction. Thus, the total projection of the spin along z, 
i.e., = X^iLi 1 is a conserved quantity in the absence of coupling to the substrate. 
In other words, the number of magnons 6^6^ is a conserved quantity of Hamiltonian 
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in (20l. However, in tlie presence of Hj^ [Eq. (24)], is no longer conserved (i.e. 
the magnons bi can "leak out" from the chain), and therefore total magnetization 
fluctuations also become damped out by the coupling to the substrate. 

The above treatment is valid provided the coupling to the substrate is weak and 
can be treated perturbatively. At stronger coupling, we have to worry about the 
possibility that the magnetic moments of the impurities are Kondo-screened by the 
substrate electrons. To deal with this situation, we shall rely on the non-perturbative 
treatment discussed in the [Appendix A[ As shown there, in the strong-coupling limit 
where Kondo correlations are important, it is possible to map the model introduced in 
Section [2] to a dissipative quantum Ising model. When coarse-grained over distances 
^ tto, this model can be described by a dissipative cj)* field theory, whose partition 
function is Z = J [d(j)] e~^^'^\ where the action S [(j)] is given by the expression |11| : 

= |<^(9,c.„)l' + | f dxdr<l>*{x,T), (39) 

being </)(x,r) = {Lhj3y^''^Y.q,uj„'i^{(l^^n) e'f'?^^""^) (/3 = 1/T) is the coarse-grained 
magnetization along the easy axis, i.e. {(j){x — Xi)) = R^^ J2fL'^R/2 i^i+j) ^ where 
1 ^ i? <C iV. This result suggests that the quantum critical point (QPC) separating 
the paramagnetic and ferromagnetic ground states (cf. Fig|4]) can be described by this 
field theory, and therefore the critical exponents around the QCP will correspond to 
those of the Wilson-Fisher fixed-point of the renormalization group (see e.g. Ref. |56| . 
for a survey). However, in deriving the map to the dissipative quantum Ising model, 

we have neglected the Heisenberg term ~ ^ J2i [^i'S''^^ + H.c] from Eq. ( 1 1. If this 
term is included, we cannot rule out that the actual values of the critical exponents 
near the QCP turn out to be different from those of the Wilson-Fisher universality 
class. Addressing this question in detail require a numerical investigation of the full 
model, which is beyond the scope of the present work. 

Nevertheless, the field theory of Eq. (39 1, is capable of describing the phases of the 
system. Focusing on the ferromagnetic phase where (^(x,t)) ^ 0, which corresponds 
to r < 0, it can be seen that the magnon excitations are indeed described by the 
propagator: 

(0*(g,w„)0((7,t^„)) = , 1 , n 2^ F ' ("^^^ 



where Eg — 2\r \ ^ Jfj — is the energy gap. This magnon propagator indicates that 
the magnetic excitations of the chain in the strong coupling regime are overdamped by 
their coupling to the metallic substrate. Physically, this makes sense as the impurity 
spins are screened by the metal electrons and thus their magnetic moments undergo 
strong damping by the collective nature of their screening clouds. This regime must be 
contrasted with the weak-coupling regime, where we found the damping to be weak. 
Thus, a physical picture emerges that allows us to distinguish the two mean-field 
solutions that we briefly described at the end of Section [3] (cf. Fig. |3|: In the weak 
coupling regime, the strong Weiss field (i.e. h^g 3> T^) due to the FM coupling of 
the impurities along the chain does not allow for the Kondo screening of the local 
moments of the atoms in the chain. As a consequence, the magnetic excitations of the 
chain are only weakly affected by their coupling to the substrate, and thus acquire a 
small line-width . On the other hand, in the strong coupling regime, the Weiss field is 
much weaker (i.e. h^s <^ Tj^) and the magnetic moments are effectively screened by 
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Figure 5. Spin-response Im Xg i^) ^ ^ function of q and u). This quantity is 
proportional to the SPEELS d P/dfldu! signal and provides information about 
the magnon dispersion relation (red line in the bottom q — ui plane) and the FM 
Ising gap. The broadening of the magnon excitation is originated in the (Kondo) 
coupling to the metallic substrate Eq. ([2|. 

the metal electrons, leading to overdamped magnetic excitations. In the next Section, 
we shall study some experimental consequences of our findings. 



5. Relation to SPEELS experiments 



In this Section, we focus on the efi^ects of the dissipative environment on observable 
quantities. The magnetic properties of low dimensional spin systems deposited 
on metals have been studied with a variety of experimental techniques, like X-ray 
magnetic circular dichroism (XMCD), used to measure directly the magnetization 
of FM ID spin chains j^l], or local probes, like spin polarized STM, which can 
provide information on the magnon dispersion relation [57J. Here we focus on the 
spin-polarized electron energy-loss spectroscopy (SPEELS) experiment |58H60| . which 
provides direct information of the magnon dispersion relation. 

The SPEELS cross section (PP/dftdut corresponds to the fractional number of 
electrons which emerge onto the solid angle dfl after being scattered by a magnetic 
excitation ui (q) in the energy range duj. According to the theory of the SPEELS 
experiment [58 , d^P/dVlduj is related to the spin response function by 



cx - Im Xq (w) , (41) 



dQ,djjj 

where (w) is the retarded spin correlation function obtained by analytic 
continuation from Xq (t^) = [Xq (^")] ^,^^^+^0+'^^*^^ 

X± (c.„) = j d^^j dr e'l-^e-^""- (5" (R, r) S+ (0, 0)) . (42) 
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Using the HP representation Eqs. ( 17 1-(|19[), we can express 
xf « 25" {bg (w„) bg (ujn)) , 

ihuJn - a\uJn\ - Eg ' 



(43) 
(44) 



where we have used Eqs. (36) and (37), and the fact that in the ID geometry q g. 
Introducing this result into Eq. ( 41 1 and performing the analytic continuation to real 
frequencies, we obtain: 



Im xf i^) = - 



(/iw - EgY + a^uj^ 



(45) 



In Fig. ^ we show the spin response function — Im {^) as a function of q and uj. 
The curves — Im x^ {^) show resonances centered at the magnon frequencies lu ~ Eg 
(red line in the bottom q — uj plane), which are broadened by the effect of a, the 
dimensionless coupling to the metal. 



6. Conclusions and Outlook 



We have studied the phase diagram and the excitation spectrum of a magnetic chain 
of atoms or molecules with easy-axis ferromagnetic interactions. We focused on a 
simple model of S* = ^ magnetic impurities displaying a Kondo-exchange interaction 
with the substrate. In the Ising limit, which provides a good approximation to the 
ground state of an easy-axis ferromagnet, we obtained the phase diagram using a 
mean-field theory approach. We find that this system exhibits two possible phases 
at zero temperature: a paramagnetic phase where the impurity spins are screened 
by the substrate, and a gapped ferromagnetic phase whose excitations are damped 
by the magnetic (Kondo) exchange with the metallic substrate. Using a bosonic 
representation for the spins as well as various types of mathematical mappings between 
models, we have also investigated the excitation spectrum in the ferromagnetic phase, 
which may be accessible through probes like spin-polarized electron energy loss 
(SPEELS). Although we have focused mainly on an easy-axis ferromagnetic (FM) 
chain, it is worth describing how our results can be modified in the case of the 
anti-ferromagnetic (AFM) chain corresponding to the parameter regime A < — 1 
of the XXZ chain described by Eq. Q. An experimental motivation to extend 
our results to this regime can be found, for instance, in Ref. [24 . Relying on 
the local bath approximation (LB A) and in the Ising limit, it is possible to apply 
the transformation S"? {—lySf, which maps the AFM to onto FM Ising chain. 
In presence of the Kondo exchange, this transformation must be supplemented by 
another rotation which takes (Ri) (—1)' (R-i) and (Ri) O s** (Ri) provided 
that (—1)' = —1. For A < —1, the excitations are separated from the ground 

state by a gap ^ ^ (Jh)^ ^ {'^h)^ ■ Using the Holstein-Primakoff representation |51] . 
we can compute the spectrum of an antiferromagnetic XXZ chain, which yields 

^AFM _ 2S\J {Jfj)^ ~ [J-^ cosqao)^ . The results of Sections |4j and js] carry over to 
the AFM (A < -1) case provided we replace Eg by E^™. A feature of the AFM 
regime that is also worth noticing is that bound states of magnons are absent from 
the spectrum [251 ES] . 

Finally, let us discuss our outlook for future research directions. One such 
direction is the study of the phases of chains of impurities with higher (i.e. S > ^) 
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spin. In this case, it is also possible to apply the same mean-field theory employed in 
Section [sj However, in order to obtain the impurity magnetization mimp {hef{,T) = 
(S'^) we must in general resort to an impurity solver such as Wilson's numerical 
renormalization group |16l EU| . This is specially true for S" > which will will 
undergo a several stage Kondo eff'ect [121 ES] and possibly underscreening jM] . Another 
interesting direction will be to asses the accuracy of the local bath approximation 
employed throughout. This remains quite a challenge, which will most likely require 
to rely on numerical methods, such as quantum Monte Carlo. Indeed, it is possible 
that this approximation only provides an rough picture of the paramagnetic phase, 
which must be corrected by the inclusion of the coupling between the different baths 
that screen the magnetic impurities individually. However, we believe that such effects 
will be less important in the ferromagnetic phase (specially for large Jfj), due to the 
existence of a gap that separates the ferromagnetic ground state from the excitations 
and therefore protects the ground state from perturbations. 

Appendix A. Bosonization of the fermionic chains and mapping to the 
dissipative quantum Ising model 

In this Appendix, we implement the Abelian bosonization approach to the semi-infinite 
fermionic ID chains along the y direction in Fig. (|2| pil This procedure is 

standard and has been successfully applied to describe the low-energy properties 
of the single Kondo-impurity. We refer the reader to the standard bibliography 
|3H l43l mi f62M64| . In the present case, the LBA greatly simplifies the complexity 
of the problem and enables a straightforward generalization to the case of many- 
impurities, each one coupled to an independent fermionic bath. At low temperature, 
in the bosonic representation [51] the Hamiltonians (|3| and ([2| become 



Hk = 



E 

i/={c,s} 

26.. 



47r 



St 



) 

V0f^, (0) 



(A.l) 



2 



S 



27r&o 



(0) 



5*; 



giy205(O) 



(A.2) 



where (jif^ {y) , (pf^ (y) are bosonic chiral fields which obey the commutation relations 
[(/>f^j, (y) ) 4'f,r] (.y')] = iTTsign (y — y') 5i^j5v,r^, and which are related to charge and 
spin density-fluctuations in the ID fermionic chains through the relations pi [y] = 
-iV(/)5 {y) and Si (y) = -^V0f, (y), respectively [HI EH- In Eq. ( [a1] vf is the 
Fermi velocity, and in Eq. (A.2 1 Sg = tan~^ (irpoJ^ba/A) is the scattering phase-shift 



associated with the potential ,1^8^/2, po = {2'kvf) the conduction electron density 
of states at the Fermi energy, and fog the lattice parameter in the fermionic chain. 
For simplicity we assume these parameters to be identical for all chains. We then 
introduce the (Emery-Kivelson) unitary transformation |43| 



U — exp 



(0) 



(A.3) 
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under which the bosonic field V(/)^^ (y) and the spin operator transform as 

Z^t (y) U = (y) + S (y) 27rjS!] , (A.4) 

U^S+U =5+e'-^*-(°). (A.5) 

Upon this transformation, the total Hamiltonian of the spin-chain coupled to the 
metallic bath, H = Ho+ Hk + Hp, transforms as H = W HU = Hq + Hk + Hp, with 



Hn = - 



H 



K 




V0f,(O) 



(A.6) 



^/27r 

g-»(V2-7)C,(0) 



H.c. 



'0 



Hp — Hj^ 



(A.7) 
(A.8) 

where we have defined Sg = Sg — 777/2 \/2. Note that the local bath approximation 
Eq. (|3| is crucial to implement bosonization along the chains, and to put 
these ideas on a clear mathematical framework. It is also interesting to note 
that in the transformed representation, the quantum dynamics of the bath 
[represented by the chiral field (0)] appears explicitly in the Heisenberg term 

(^5+S'r^ie*''['^^=(°)-"^'"+i.='"^] +H.C.) [15,. Physically, this means that the 

Heisenberg interaction is now "dressed" by the spin-density fluctuations of the electron 
gas. In the case of an easy-axis spin chain in Ising limit Jfj S> J^, the effect of this 
term is negligible and can be ignored in a first approximation. 

We now exploit the freedom to choose 7 and set 7 = ^/2. In this case, the 
Hamiltonian reads 



H — Hn + Hi, 



Hp 



E 



JpfSj s. 



2~6s '^<t>t (0) 



71" PO 



'2ti 



2tt ' 



Hf 



(A.9) 



This Hamiltonian corresponds to the ID dissipative quantum Ising model (DQIM), 
where now the Kondo Hamiltonian H^ is equivalent to the spin-boson model with 
Ohmic dissipation [U |5J [T31 [S3] with 5s related to the dissipative parameter a in 

the context of macroscopic quantum coherence through a = [25s/t^ , and with 
the in-plane Kondo interaction playing the role of a magnetic field along the a;— axis 
hx = — Jj^/27r. At r = 0, the ID DQIM is known to display a paramagnetic to 
ferromagnetic quantum phase transition which is in the universality class of quantum 
dissipative systems and whose dynamical critical exponent is z = 2 [11^. Note 
that this is very different with the case of the ID quantum Ising chain, which is 
in the universality class of the 2D classical Ising model and where z = 1. The 
critical properties of this model near the quantum phase transition have been studied 
in the context of antiferromagnetic instabilities of Fermi liquids |lll I56| using the 
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framework of the Hertz-Moriya-Millis theory |65fl67j This theory describes critical 
quantum fluctuations of the order parameter around the Gaussian fixed point of the 
theory. The predicted value of the critical dynamical exponent z = 2 and has been 
confirmed numerically with Monte Carlo simulations I13|. Physically speaking, 
such a dynamical exponent implies that the effective dimensionality of the ID spin 
chain coupled to the metallic bath is des = I + z = 3, and therefore fluctuations 
of the order parameter are expected to be much less important than the case of 
the non-dissipative classical or quantum Ising model. This fact supports our MF 
approximation, which allows to extend these results to T > 0. 
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